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■ ABSTRACT 

, Aims. We examine whether dust mass can trace the total or molecular gas mass in late-type Virgo cluster galaxies, and how the 

environment affects the dust-to-gas ratio and the molecular fraction. 
Q ' Methods. Using the far-infrared emission, as observed by the Herschel Virgo Cluster Survey (HeViCS), and the integrated HI and 

5_( ' CO brightness, we infer the dust and total gas mass for a magnitude limited sample of 35 metal rich spiral galaxies. Environmental 

, disturbances on each galaxy are considered by means of the HI deficiency parameter 

■ Results. The CO flux correlates tightly and linearly with far-infrared fluxes observed by Herschel, especially with the emission at 
160, 250 and 350//m. Molecules in these galaxies are more closely related to cold dust rather than to dust heated by star formation or 
to optical/NIR brightness. We show that dust mass establishes a stronger correlation with the total gas mass than with the atomic or 
molecular component alone. The correlation is non-linear since lower mass galaxies have a lower dust-to-gas ratio. The dust-to-gas 

^ . ratio increases as the HI deficiency increases, but in highly HI deficient galaxies it stays constant. Dust is in fact less affected than 

OO ' atomic gas by weak cluster interactions, which remove most of the HI gas from outer and high latitudes regions. Highly disturbed 

04 ' galaxies, in a dense cluster environment, can instead loose a considerable fraction of gas and dust from the inner regions of the disk 

, keeping constant the dust-to-gas ratio. There is evidence that the molecular phase is also quenched. This quencing becomes evident 

. by considering the molecular gas mass per unit stellar mass. Its amplitude, if confirmed by future studies, highlights that molecules 

' are missing in Virgo HI deficient spirals, but to a somewhat lesser extent than dust. 

' Key words. Galaxies: clusters - ISM: Molecules - infrared: dust,galaxies - individual: Virgo 
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1. Introduction relation between dust and molecules it is impoitant to have 

a good inventory of their masses and distribution in a vari- 

Dust and gas in galaxies are essential ingredients for star for- g^j^^y typ^^ environments. Intense interstellar radia- 

^ • mation. Stars are born because of the coohng and fragmenta- ^^^^ ^gj^^ galaxy-galaxy interactions or gas stripping as the 

. . tion of molecular gas, which in today's galaxies forms from g^j^^y ^^^^^ ^^rough inter-cluster gas can quench the growth 

atomic gas pnmaril^b^talytic reactions on the surface of molecular clouds through galaxy starvation or strangulation 

dust grains ( |Gould & Salpeteij|1963|) . Dust and atomic gas are, (Larson et al. 1980; Bekki et al. 2002; Boselli & Gavazzi 200^ 

however, only part of the picture because the fractional abun- iFumaealli & Gavazziii2008i; .Fumagalli et al. 2009). But there is 

dance of molecules depends on a delicate balance between for- ^ unanimous consent on how environmental conditions affect 

mation and dissociation processes (|Hollenbach et al. ||I971) and j^e molecular component: this is more bounded to the galaxy 

on the self-shielding ability of molecular species from ultravio- potential well, being more c onfined insi de the disk, and less dis- 



let radiation (Draine & Bertoldi 1996). The presence of promi- turbe d bv the environment (IStark et al.l flglet ISolomon & Sagil 

nent arms and high interstellar pre ssure i n the disk favour high llQgsLlKennev & Youngll 19891) 
molecular abundances ( Elmegreeiil l 1 9931) . In order to study the 

Measuring the total molecular content of a galaxy is a rather 

* Herschel is an ESA space observatory with science instruments difficult task because molecular hydrogen (H2), the most abun- 

provided by European-led Principal Investigator consortia and with im- dant molecular species, is hard to excite, lacks a permanent elec- 

portant participation from NASA. trie dipole moment and must radiate through a slow quadrupole 
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transition. The ability of molecular hydrogen to form stars, the so 
called star formation efficiency, can also vary and hence the star 
formation rate cannot be used a s an i ndicator of mo lecular mass 
(lYounglll999t iDi Matteo et"aLll2007t lMurravll201 II) . Molecular 
surveys rely on the brightness of CO rotational lines and on CO- 
to-H2 conversion factors, Xco, to derive the molecular hydro- 
gen mass from the CO luminosity. In general this depends on 
the CO-to-H2 abundance ratio, which is a function of metallic- 
ity and on the self-shielding ability of the two molecules, and 
on the excitation state of the CO molecules. The CO molecule 
forms almost exclusively in the gas phase and has a lower ca- 
pability for self-shielding then H2 due to the low abundance 
of carbon relative to hydrogen (e.g. iGlover & Mac Low 2011 ). 
However, there are still discrepant results on how Xgo varies, for 
example as a function of the gas metallicity ('Wilson" 1 995'; 'Israel' 
[l997b; Bolatto et al. 2008; Leroy et al. 2011; Shetty et al. 2011). 
Infrared emission has been used sometime t o find possible vari- 
ations of the CO-to-H2 conversion factor (iKnapp et alj T987t 



Israelll997allbl:rLerov et a .l2007HGratier et al.l2010llLerov et alJ 



2011t|Magnni et alj|201 1|) and estimate the molecular content. 

The total dust mass of a galaxy can only be determined if 
dust emission is traced throughout using multi-wavelength ob- 
servations from the infrared (hereafter IR) to the sub-millimeter. 
Hot dust, in the proximity of star-forming regions, is similar to 
the tip of an iceberg: prominent but not representative of the to- 
tal dust mass. The Hersche l Space Observatory ( Pilb ratt et al.l 
I2OIOI; iPoplitsch et all 1201 Ot iGriffin et alj 1201 Ol) . with its two 
photometric instruments at far infrared wavelengths (PACS, the 
Photodetector Array Camera Spectrometer with 70, 100, and 
160 yum bands and SPIRE, the Spectral and Photometric Imaging 
Receiver with 250,350 and 500 /vm bands) can observe both 
sides of the peak (100 - 200jum) in the spectral energy distri- 
bution of dus t in galaxies. The Herschel Virgo Cluster Survey 
dPavies et al.P 2010. hereafter HeViCS ), an Herschel Open Time 
Key Project, has mapped an area of the nearby Virgo cluster 
of galaxies, with PACS and SPIRE to investigate the dust con- 
tent of galaxies of different morphological types and whether 
intra- cluster d ust is present. Recent studies based on Herschel 
data dCortese et al. 2010, 2012) have already pointed out that in 
a dense cluster environment dust can be removed from galaxy 
disks. 

For a sample of bright spiral galaxies in the Virgo cluster, 
which have been targets of CO observations, we shall examine 
the correlation between the dust emission and the CO bright- 
ness, integrated over the disk, and between dust and gas masses. 
An important question to address is, for example, whether HI 
deficient galaxies are also deficient in their global dust and 
molecular content and ho w deficiency affects galaxy evolution 
(iBoselli & Gavazzil l2006l) . The large number of Virgo spiral 
galaxies, forming stars actively, will let us investigate whether 
dust mass establishes a better correlation with the molecular gas 
mass or with the sum of the atomic and molecular gas mass and 
to pin down possible environmental effects on the dust-to-gas 
ratio. A large sample of field galaxies would be an ideal refer- 
ence sample to have, in order to distinguish environmental ef- 
fects from properties related to the galaxy quiescent evolution. 
However, the presence of newly acquired cluster members in a 
young cluster, such as the Virgo Cluster, offers the possibility 
of using the HI deficiency as indicator of environmental distur- 
bances. By defining the HI deficiency as the difference between 
the observed HI mass and that expected in an isolated galaxy, 
[Oiovanelli & Havnes ( 1985) have in fact shown that the HI defi- 
ciency increses towards the cluster center This supports the idea 
that the HI deficiency can trace environmental effects and that 



galaxies in the same cluster which are not HI deficient can be 
considered similar to isolated galaxies. 

Another open question is whether dust-to-gas ratio can be 
used as metallicity indicator If the production of oxygen and car- 
bon are the same and the production, growth and destruction of 
solid grains does not change from a galaxy to another, we expect 
the fraction of condensed elements in grains to be proportional 
to those in the gas phase. Radial analysis of i ndividual galax- 
ies has shown that this is indeed the case (Ma grini et aLlbofll ; 
iBendo et 301201 Obh and models have been proposed to predict 
how the dust abundance varies as galaxies evolve inmetallicity 
(Lisenfeld & Fesaral [HH iHirashita et a"lll2002b . iDraine et alJ 
(2007) found that the dust-to-gas ratio in the central regions does 
not depend on the galaxy morphological type for galaxies of type 
Sa to Sd. But a linear scaling between the metallicity and the 
dust-to-gas ratio is not always found (Boselli et al. 20 0^, es- 
pecially when metal abundan ces are low ( Muiioz- Mateos et al.l 
12009) . However, as shown bv lGalametz et al. (201 L) , this might 
result from an inaccurate determination of dust mass due to the 
lack of data coverage at long infrared wavelengths. Using the 
wide wavelength coverage of Herschel we investigate whether 
the dust-to-gas ratio can be used as metallicity indicator for 
metal rich spirals in Virgo. 

We describe the galaxy sample and the database in Section 2, 
the relation of the CO brightness and gas masses with the far- 
infrared emission and dust masses in Section 3. An analysis 
of possible environmental effects on the dust and gas content, 
atomic and molecular, is given in Section 4. In Section 5 we 
summarize the main results and conclude. 



2. The sample and the data 

Our sample is a magnitude-limited sample of galaxies in the 
Virgo cluster area mapped by HeViCS. The HeViCS covers 
=a 84 sq. deg. of th e Virgo cluster, not the whole cluster area 
dPavies et alj 1201 Ol) . and the full-depth region observed with 
both PACS and SPIRE is hmited to 55 sq. deg., due to offsets 
between scans. Galaxies in the Virgo Cluster Catalogue (VCC) 
(Binggeli et al. 1985) are included in our sample only if they are 
of morphological type Sab to Sm, are located in the HeViCS 
fields, their distance is 17 < D < 32 Mpc, and their total mag- 
nitude in the B-syst em, as given in the Third Re ference Catalog 
of Bright Galaxies (Ide Vaucouleurs et al.l ll995). is BT< 13.04. 
This limiting magnitude ensures that all galaxies in our sample 
have been observed and detected at 21 -cm and in the CO J=l-0 
line. Very few Virgo galaxies fainter than BT= 13.04 have been 
detected in the CO J=l-0 line. In the HeViCS area, and for the 
selected range of morphological type, our sample is 100% com- 
plete down to B-mag 13.04. With respect to the 250 jum flux our 
sample is 100% complete down to F250 =5 Jy. The significance 
of correlations between different galaxy properties are assessed 
using the Pearson linear correlation coefficient, r. 

The complete list of 35 galaxies selected in our sample is 
given in Table 1 . In Column 1 we give the sample code, which 
indicates the type of CO observations available for estimating 
the molecular mass. "M" refers to the main sample, which con- 
sists of all galaxies in HeViCS that have been fully mapped in 
the CO J=l-0 line. "S" refers to the secondary sample: galax- 
ies in HeViCS which have been detected but not fully mapped 
in CO. The galaxy name, right ascension and declination are 
given in Columns 2, 3, and 4 respectively. The distance to the 
galaxy in Mpc, the galaxy morphological type and stellar mass 
are given in Columns 5, 6, and 7, respectively. Distances are de- 
termined from the average redshift of the individual cluster sub- 
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groups (iGavazzi et al ] ll999h and we adopts the morphological 
classification scheme of the VCC. Nebular oxygen abundances, 
12+log(0/H), from the Sloan Sky Digital Survey are shown in 
Column 8. In parenthesis the metallicity inferre d from the stel- 
lar ma ss using the mass-metallicity relation of iTre monti et alJ 
( 2004 ) and adopting the functional form fit of lKewle v & Elliso^ 
poo s') . The atomic and molecular hydrogen mass are listed in 
Columns 9, 10-11, respectively. The two values of the molecu- 
lar mass, and , have been derived from the integrated 
CO brightness using two different values of the CO-to-H2 con- 
version factor: a uniform constant value of Xco or a metallicity 
dependent CO-to-H2 conversion factor (as given in more details 
in Section 2. 1). In Column 12 we give the relevant reference for 
the integrated CO line intensity. In the last column we indicate 
the number of positions where CO has been detected/observed 
or whether the galaxy has been fully mapped. If more than one 
map is available we take the mean CO brightness and list the 
references joined by a plus sign. 

We derive the stellar mass of each galaxy, M ,, using the B-K 
color and the K-band luminosity, as described in Ma nnucci et alj 
(l200l . In order to have a uniform set of photometric data, 
we use the GOLDMINE databasf] (Gavazzi et al. 2003) for th e 
Ha, B-, H- and K-band magnitudes ( Gavazzi & Boselli 
and for the disk geometrical parameters (Nilson 1973). We de- 
rive absolute luminosities after correcting fo r internal extinc- 
tion a ccording to the prescription given by Gavazzi & Bosellil 
(119961) . Corrections for foreground galactic extinction are small 
(Av ^ 0.07). 

As a metallicity ind i cator we use the O/H abundance ob- 
tained bv lTremonti et all ( |2004 from the SDSS data. Metallicity 
has been estimated by fitting at the same time all prominent 
emission lines ([O II], H^, [O III1 Hg, \N 111, fS 111) with a n 
integrated galaxy spectrum model (ICharlot & L onghett? l2001h . 
Not all galaxies in our sample have SDSS spectra with measur- 
able emission lines or appropriate line ratios needed to derive 
O/H abundances. This is because SDSS narrow slits are gen- 
erally positioned on the central high surface brightness regions 
and the presence of some nuclear activity or of a bulge, pre- 
vents the authors from probing emission line HII regions in the 
disk. Galaxies with very weak emission lines or with anoma- 
lous line ratios because of nuclear activity have been discarded 
by Tremonti et al. (2004). Thanks to the analysis of large galaxy 
samples from the SDSS, the mass-metallicity relation predicts 
that metallicity scales with stellar mass until the stellar mass is 
below ~ 3 X 10'" M0, then it flattens out. Our sample is repre- 
sentative of the mass-metallicity relation: the dispersion of the 
original dataset around the functional form fi t is ^ 0.15 dex 
(ITremonti et all 120041: iKewlev & EUisonl l2008h . which is con- 
sistent with that of our sample. In this paper we shall use the 
O/H abundance derived from the mass-metallicity relation for 
all galaxies with known K-band luminosity. 

Values of atomic hydrogen mass are taken from the VLA 
Imaging of Virgo Spirals in Atomic Gas database (VIVA, 
IChung et alj |2009a). where possible. We correct the HI masses 
to reflect the galaxy distances used in the present paper if they 
differ from t he value of 17 Mpc adopted in the VIVA survey. 
IChungetalj|2009 a) have found an excellent agreement between 
single dish fluxes and VIVA fluxes, with no indication that VIVA 
data is missing any very extended flux. If the galaxy is not in 
VIVA the GOLDMINE value of the HI mass is used. We wiU 
use Mgas to indicate the total hydrogen mass of a galaxy i.e. the 
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sum of the atomic (M^/) and molecular hydrogen mass (M/^J 
corrected for helium abundance. 



2.1. Integrated CO J=1-0 brightness and molecular mass 
estimates 

Extensive searches for CO emission in Virgo galaxies have 
been carried out by a number of authors using single dish mm- 
telescopes. Only 18 galaxies in the HeViCS area have been 
mapped (M-sample, plotted with filled circles). Some galaxies 
have been mapped by more than one team using different tele- 
scopes. In this case the total CO fluxes are in good agreement and 
we use the arithmetic mean of the two values. Uncertainties from 
the original papers have been adopted. We include in the M- 
sample NGC4522 which has been mapped in the CO J=2-l line 
and observed at several positions in the CO J=l-0 to determine 
the 2-1/1-0 line ratio. We do not include NGC4299 whose CO 
map has a low signal to noise and allows a reliable determina- 
tion of the total CO brightness only within the Ha emission map 
boundary. There are 17 galaxies with BT< 13.04 which have 
been detected in the CO J=l-0 line with one or several pointed 
observations across the optical disk (S-sample, plotted with open 
circles). All but one of these have CO J=l-0 line emission above 
3-0- in at least one beam. Only for NGC4206 the CO J=l-0 line 
has been marginally detected, at the 2-cr level. When the CO line 
brightness has been measured in more than one position we rely 
on the total CO flux estimated in the original papers. In these 
cases the total flux is determined via exponential or gaussian fits 
to the radial distribution (Stark et al. 1986; Young et al. 1995). If 
a galaxy has been detected at just one position we determine its 
total CO flux assuming an exponential radial distribution, cen- 
trally peaked, with scalelength A = 0.2/?25- 

The value A - 0.2 + 0.1/? 25 is the average radial scalelength 
of the CO flux distribution inferred for the M-sample (mapped 
galaxies). This is consistent with the exponential scalelengths 
measured for the CO emission in nearby galaxies mapped in 
the Heterodyne R eceiver Array CO Line Extragalactic Survey 
dLerov et al.ll2009l) . We convolve the exponential function, cor- 
rected for disk inclination, with the telescope beam, and compare 
this to the observed flux in order to determine the central bright- 
ness and hence the total CO flux. Notice that S-sample galaxies, 
as well as mapped Virgo galaxies, do not have enhanced CO 
emission in the center and are very well r epresented by expo- 
nential functio ns, as shown by Kuno et al.l (2007). For HeViCS 
galaxies in the iKuno et alJ (12007 ) sample the fraction of molec- 
ular gas in the inner regions, f,>,, is only 0.14 on average, for 
both Seyfert and non-Seyfert galaxies. The dispersion between 
the total CO flux as measured from the map, and the value de- 
rived from the exponential fit with scalelength A = 0.2/? 25 is 0. 15 
in log units. This is the error which we add in quadrature to the 
flux measurement error when computing the uncertainties in the 
total CO brightness of galaxies in the S-sample. 

For 6 galaxies in Table 1 we use newly acquired data from 
observations with the IRAM 30-m telescope (Pappalardo et al. 
2012). Here we give a brief description of the observations. We 
searched for CO J=l-0 and J=2-l lines in NGC4189, NGC4294, 
NGC4298, NGC4299, NGC4351, NGC4388, during June 2010 
and April 201 1. The FWHM beam is 22 arcsec at 1 15 GHz, the 
frequency of the J=l-0 line used in this paper We observed the 
sources in position-switching mode, using the EMIR receiver 
combination E0/E2 and the VESPA and WILMA backend sys- 
tem with a bandwidth of 480 MHz and 4 GHz, respectively. The 
on-the-fly mapping mode was used for mapped galaxies. Non- 
mapped galaxies were sampled along the major axis at 22 arcsec 
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Table 1. Galaxies in the M and S samples and related quantities. 





ID 


RA 


DEC 


D 


T 


M. 


O/H 


log M„, 


log MJ,^ 
Mo 


logM^^ 


Ref. 


Notes 






deg 


deg 


Mpc 




Mo 


dex 


Mo 


Mo ' 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


M 


NGC4189 


183.447 


13.425 


32 


Sc 


12.69 


9.1(9.1) 


9.30" 


9.26 


9.37 


1 (2) 


map 


M 


NGC 4192 


183.451 


14.900 


17 


Sb 


10.73 


(••) (9.1) 


9.63° 


9.05 


8.90 


3 


map 


M 


NGC 4212 


183.914 


13.901 


17 


Sc 


11.82 


9.1(9.1) 


8.91 


9.05 


8.96 


3 


map 


M 


NGC 4254 


184.707 


14.416 


17 


Sc 


10.45 


9.2(9.1) 


9.65° 


9.88 


9.60 


4+3 


map 


M 


NGC 4298 


185.387 


14.606 


17 


Sc 


11.95 


9.1(9.1) 


8.69° 


9.02 


8.96 


4+1 


map 


M 


NGC 4302 


185.427 


14.598 


17 


Sc 


12.31 


(••) (9.1) 


9.17° 


9.09 


8.93 


4 


map 


M 


NGC 4303 


185.479 


4.474 


17 


Sc 


10.30 


9.2(9.1) 


9.68 


9.76 


9.49 


4+3 


map 


M 


NGC 4321 


185.729 


15.823 


17 


Sc 


10.02 


9.2(9.1) 


9.46° 


9.79 


9.39 


4+3 


map 


M 


NGC 4388 


186.445 


12.662 


17 


Sab 


11.87 


8.8(9.1) 


8.57° 


8.85 


8.90 


1 (5,2) 


map 


M 


NGC 4402 


186.531 


13.113 


17 


Sc 


12.64 


9.0(9.1) 


8.57° 


9.16 


9.10 


4+3 


map 


M 


NGC 4438 


186.940 


13.009 


17 


Sb 


11.12 


(••) (9.1) 


8.68 


8.83 


8.56 


4 


map 


M 


NGC 4501 


187.997 


14.421 


17 


Sbc 


10.50 


8.8(9.1) 


9.22° 


9.75 


9.90 


4+3 


map 


M 


NGC 4522 


188.415 


9.175 


17 


Sbc 


12.97 


9.0(9.0) 


8.53° 


8.48 


8.51 


5 


map 


M 


NGC 4535 


188.585 


8.198 


17 


Sc 


10.73 


(••) (9.1) 


9.52° 


9.50 


9.25 


4+3 


map 


M 


NGC 4567 


189.136 


11.258 


17 


Sc 


11.91 


(..) (9.1) 


8.97° 


9.05 


8.97 


4 


map 


M 


NGC 4568 


189.143 


11.239 


17 


Sc 


11.22 


9.2(9.1) 


9.18° 


9.35 


9.03 


4 


map 


M 


NGC 4569 


189.208 


13.163 


17 


Sab 


10.08 


(..) (9.1) 


8.79° 


9.46 


9.16 


4+3 


map 


M 


NGC 4579 


189.431 


11.818 


17 


Sab 


10.52 


(••) (9.1) 


8.75° 


9.18 


9.04 


4+3 


map 


S 


NGC 4152 


182.656 


16.033 


32 


Sc 


12.77 


9.3(9.1) 


9.73 


9.20 


9.13 


6 


1/1 


S 


NGC 4206 


183.820 


13.024 


17 


Sc 


13.00 


(..) (9.0) 


9.38 


8.47* 


8.50* 


6 


2/5 


S 


NGC 4216 


183.977 


13.149 


17 


Sb 


10.84 


(..) (9.1) 


9.25° 


9.10 


9.03 


2 


4/9 


S 


NGC 4237 


184.298 


15.324 


17 


Sc 


12.66 


(..) (9.1) 


8.32 


9.15 


9.08 


6 


1/5 


S 


NGC 4273 


184.983 


5.343 


32 


Sc 


12.49 


9.1(9.1) 


9.54 


9.32 


9.25 


2 


2/3 


S 


NGC 4294 


185.324 


11.511 


17 


Sc 


12.68 


(••) (9.0) 


9.21° 


7.85 


7.88 


1 (2) 


2/2 


S 


NGC 4299 


185.420 


11.503 


17 


Scd 


12.99 


8.8(8.6) 


9.04° 


8.08 


8.53 


1 (2) 


map 


S 


NGC 4307 


185.523 


9.044 


23 


Sbc 


12.78 


(..) (9.1) 


8.15 


8.97 


8.90 


7 


3/3 


S 


NGC 4312 


185.630 


15.538 


17 


Sab 


12.46 


9.1(9.0) 


8.08 


8.87 


8.91 


2 


1/3 


S 


NGC 4313 


185.660 


11.801 


17 


Sab 


12.50 


(..) (9.1) 


8.02 


8.82 


8.75 


7 


1/3 


S 


NGC 4351 


186.006 


12.205 


17 


Sc 


12.99 


9.0(8.9) 


8.48° 


7.86 


8.00 


1(5) 


1/1 


s 


NGC 4380 


186.342 


10.017 


23 


Sab 


12.43 


(..) (9.1) 


8.37° 


8.76 


8.69 


8 


1/1 


s 


NGC 4413 


186.634 


12.611 


17 


Sbc 


12.94 


9.1(9.0) 


8.29 


8.43 


8.46 


9 


1/1 


s 


VCC 939 


186.697 


8.885 


23 


Sc 


13.06 


(..) (9.0) 


9.34 


8.31 


8.35 


8 


1/1 


s 


NGC 4430 


186.860 


6.263 


23 


Sc 


12.76 


9.0(8.9) 


8.86 


8.76 


8.90 


9 


1/1 


s 


NGC 4519 


188.376 


8.654 


17 


Sc 


12.46 


9.1(9.0) 


9.43 


8.32 


8.35 


8 


1/1 


s 


NGC 4571 


189.235 


14.217 


17 


Sc 


12.15 


(..) (9.1) 


8.79 


8.93 


8.91 


2 


3/5 



References. (a)HI mass from the VIVA database; (ii)this mass is highly uncertain because it has been estimated using a CO flux only a factor 
two above the rms; (c) mass for a constant CO-to-H2 conversion factor; (v) H2 mass for a metallicity dependent CO-to-H2 conversion factor. 
References (Column 12). (l )new I RAM-30mt observations, C.Pappalardo et al. (2012) (in parentheses the refer ence to the older, lower sensitiv- 
ity puWished data), (2 f^ng et al.l (1 1995).(3: Kuno et al. (2007),(4; Chung et al. (2009b), (5 Vollmer etaP ( l2008l) .(6 lStark et al.l ( li986h ,(7)Bosem 
private communication (see also GOLDMINE database).(8 '.'Boseni et al.. (.2002) .(9 .BoseUi et al. {.199% 



spacing. The spectra were smoothed in velocity to 10.5 km s"' 
and the data from different backends were averaged resulting in 
an rms between 7 and 19 mK. 



oxygen abundances in our own galaxy through optical studies 
(iRudolDh et al.ll2Q06h . This can be rewritten as: 



R 

-0.8 — +Z, 

"25 



(1) 



Two values of the molecular hydrogen mass for each galaxy 
are listed in Column 10 and 11 of Table 1. In Column 10 
M^^ is derived from the total integrated CO flux using the 
given distance and a constant CO-to-Ha conversion factor, 
equal to that found in the solar neighborhood Xco - 2 x 
10^°mols. cm~^/(Kkni s~') (Strong & Mattox 1996; Da me et al.l 
l200lllAbdo et al.l2010l;IShettv et al.F201 11) . In Column 1 1 M)^^ is 
the molecular hydrogen mass for a metallicity dependent CO-to- 
H2 conversion factor computed as follows. We assume that the 
metallicity at the galaxy center, Zc, is that given by the mass- 
metallicity relation (only for VCC975 we do take the SDSS 
value because of the unknown K-band magnitude) and that there 
is a radial metallicity gradient similar to that determined by the 



where Z=12+log O/H. Using the above equation and the value 
of the CO-to-H2 conversion factor measured in the solar neigh- 
borhood (at R ^ 0.67/?25) the following radial dependence can 
be written: 



log 



^co 
2 X 102" 



0.8- 



R 

R25 



■Z, + 8.67 



(2) 



if one assumes a CO-to-H2 conversion factor which varies in- 
versely with metallicity. For Z^. = 9.1, a typical value of the cen- 
tral metallicity in our sample, the above equation predicts that 
the conversion factor is 8x10''^ cnr^/{K km s"') at the center 



and 5x10^" cm"7(K km s~') at 7? - R25- Unless stated differ- 
ently we shall use X'i„ and the corresponding masses (Column 
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11 in Table 1) throughout this paper However, due to the lack 
of data for radial metallicity gradients in cluster galaxies and to 
uncertainties of metallicity calibrators, we will quote also results 
relative to . Given the limited metallicity range of our sam- 
ple results are not muc h dependent on the CO-to-H2 factor used 
(iBolatto et aL.2008: Lerov et alJIioTTl) . 

2.2. Far-infrared observations and dust mass estimates 

For the far-infrared ( hereafter FIR) emiss ion we use the HeViCS 
dataset described in iDavies etal ] (l2012h but with all the eight 
scans combined, as described in Auld et al. (2012), since 
Herschel observations have now been completed. The Herschel 
satellites has scanned the whole area of the survey four times 
in two orthogonal directions giving a total of eight images for 
each band. The full width half maximum (FWHM) beam sizes 
are approximately 9.4 and 13.5 arcsec with pixel sizes of 3.2 and 
6.4 arcsec for the 100 and 160 //m PACS channels, respectively. 
The FWHM of the SPIRE beams are 18.1, 25.2, and 36.9 arcsec 
with pixel sizes of 6, 8, and 12 arcsec at 250, 350, and 500 /jm, 
respectively. 

As shown in iBendo et alJ (l2010ah . emission longward of 
100 fim traces the cold dust heated by evolved stars which con- 
stitutes the bulk of dust mass in spiral galaxies. We can derive 
dust masses and temperatures using the PACS and SPIRE fluxes 
after fitting the spectral energy distribution (SED) between 100 
and 500 yum. In Table 2 we show the Herschel fluxes, dust masses 
and temperatures, and the values for the SED fits relative to 
all galaxies in our sample. Dots indicate that the galaxy is out- 
side the area mapped by HeViCS in that band. The errors quoted 
in Table 2 cover uncertainties due to the choice of the aperture, 
to the calibration and to photometric random errors (Auld et al. 
2012). Calibration uncertainties are the largest errors also for 
PACS photometry, (in this case they have been assumed to be of 
order of 20%). For bright galaxies, such as those in our sample, 
the procedure and dataset described in Auld et al. (2012) gives 
fluxes, dust temp eratures an d masses which are consistent with 
those obtained bv lDavies et al. (2012) using only two orthogonal 
scans i.e. one quarter of the final exposure times. Using the com- 
plete HeViCS dataset we detect all galaxies in all 5 bands. For 
only one galaxy, NGC 4571, we cannot measure PACS fluxes 
because it lies outside the boundaries of the PACS mapped area. 

The SED fits were done for each galaxy after convolving 
the modified blackbody model with the spectral response func- 
tion for each band, thus automatically taking into account color 
corrections for the specific SED of each object. We used a 
power law dust emissivity a:,} = KoiAo/A'f, with spectral in- 
dex P - 2 and emissivity a:o - 0.192 m^ kg"' at /lo=350/im. 
These values reproduce the average emi ssivity of mo dels of 
the Milky Way dust in the FlR-submm ( lDraindl2003l) . They 
are also consistent with the dust spectral emissivity index (3 
deter mined in molecular clouds using the recent Planck data 
(e.g. iPlanck Collaboration et al.l 120111 and references therein). 
The standard minimization technique is used in order to de- 
rive the dust temperature and mass. Uncertainties in the tem- 
perature determination, considering the photometry errors, are 
about 1-2 K. Uncertainties on dust temperature and mass are 
only slightly larger for the NGC 4571 for which there is no 
PACS data. More details on SED fitting are given in lSmith et al.l 
(l2010'):'Magrini et al.l (1201 lh : lDavies et all (1201 2 *). Note that sin- 
gle black body SED models fit well S-sample galax ies as well as 
M-sa mple galaxies. All SED fits are displayed by iDavies et alj 
(12012) and by Auld et al. (2012). A two component fit is not 



well constrained by just five data points in the SED and it would 
require shorter and longer wavelength data. 

3. The correlations between dust and the molecular 
or total gas mass 

Far-infrared fluxes for galaxies in the M-sample correlate 
well with the total CO flux, as shown in Fig. [1] Pearson 
linear correlation coefficients r in the log F/r/R-log Fco 
plane are high for fluxes in all 5 Herschel bands, rang- 
ing between 0.88 and 0.93 for the whole sample, and 
between 0.94 and 0.96 for galaxies in the M-sample (the 
highest being for Fieo and F250). Slopes are of order unity 
(0.92+0.08, 1 .04+0.08, 1 .08+0.08, 1 .10+0.07, 1 .09+0.05 for 
FIR wavelengths from 100 to 500 //m), meaning that 
the correlation is close to linear The best linear fit in 
the logF25o-logFco plane gives for example logFco = 
1.08(+0.08)logF25o + 1.45(±0.08). These slopes are compatible 
with those derived considering only galaxies in the M-sample 
(1.04+0.13,1.06+0.11,1.04+0.10,1.04+0.08,1.04+0.05 for FIR 
wavelengths from 100 to 500 fim). The less tight relation 
between the CO line intensity and hot dust emission in Virgo 
spirals (e.g. IStark et al.l Il986i for IRAS data) suggests that 
molecules are more closely related to cold dust rather than to 
dust heated by star formation. 

The larger scatter and the slightly lower correlation coeffi- 
cients found when we include in the analysis galaxies in the S- 
sample, might be driven by the higher uncertainties on the total 
CO emission, because of the lack of CO maps. In order to check 
whether this is indeed the case, or whether a break in the correla- 
tion occurs at low luminosities, additional CO maps are needed 
for galaxies as dim as 13-14 mag in the B-band. To avoid the un- 
certainties on the integrated CO fluxes for non mapped galaxies 
we have checked the slope and the scatter of the FIR-CO cor- 
relation for galaxy central regions. We convolved the 250 yum 
maps from the original resolution (FWHM=18.2") to the reso- 
lution of pointed CO J=l-0 observations (which can be found in 
references to Table 1). At the optical center of galaxies we mea- 
sure the 250 ;um flux in an area as large as the FWHM of the 
telescope used to measure the CO J=l-0 line intensity (ranging 
between 45 and 100 arcsec). Figure|2]shows the resulting corre- 
lation between F250 and the CO flux in the same area, as given in 
the original papers. The straight line, of slope 1.11 +0.02, shows 
the best linear fit to these data. This is compatible with the slope 
measured when values integrated over the galaxy disk are con- 
sidered (1.04+0.08). The scatter in the correlation is not reduced 
with respect to that found for the same quantities integrated over 
the whole disk. 

Dust masses take into account temperature variations which 
individual FIR fluxes do not. The bluest galaxies have lower dust 
temperatures i.e a low F25o-to-F5oo flux ratio. In Figure |3] we 
compare the dust mass, M^,,,,, to the molecular, M//,, atomic, 
M///, and total gas mass, M^oj. The lowest scatter in the dust-to- 
molecular gas mass relation is found if M//, =MJ^ , for the whole 
sample (r = 0.81) or for the M-sample alone (r - 0.88 in this 
case). The correlation for the whole sample reads: 

logM^^ = 0.97(+0.13) X \ogMdus, + 1.9(+1.0) (3) 

A flatter, but still compatible slope is obtained if one considers 
the correlation M(/„5,-M^ : 

logM),^ = 0.77(±0.12) X \ogMdus, + 3.2(+0.9) (4) 
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Table 2. Herschel fluxes in the various bands, dust masses, temperatures and reduced for the SED fit. 



Sample 


ID 


F500 


F350 


F250 


F160 


Fioo 


log Mrf,,,, 










Jy 


Jy 


Jy 


Jy 


Jy 


Mo 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

M 

S 

S 

S 

S 

S 

S 

S 

S 

S 

S 

s 
s 
s 
s 
s 
s 
s 



NGC4189 
NGC 4192 
NGC 4212 
NGC 4254 
NGC 4298 
NGC 4302 
NGC 4303 
NGC 4321 
NGC 4388 
NGC 4402 
NGC 4438 
NGC 4501 
NGC 4522 
NGC 4535 
NGC 4567 
NGC 4568 
NGC 4569 
NGC 4579 
NGC 4152 
NGC 4206 
NGC 4216 
NGC 4237 
NGC 4273 
NGC 4294 
NGC 4299 
NGC 4307 
NGC 4312 
NGC 4313 
NGC 4351 
NGC 4380 
NGC 4413 
VCC 939 
NGC 4430 
NGC 4519 
NGC 4571 



1.02±0.08 
5.03±0.37 
1.88+0.14 
9.28+0.66 
1.87+0.14 
3.23+0.23 
8.27+0.59 
10.37+0.74 
1.38+0.10 
2.22+0.16 
1.40+0.13 
9.24+0.65 
0.59+0.05 
6.33+0.48 
1.70+0.30 
4.20+0.84 
3.27+0.24 
3.36+0.25 
0.81+0.07 
1.08+0.09 
4.01+0.29 
1.12+0.09 
1.56+0.12 
0.90+0.07 
0.49+0.05 
0.68+0.05 
0.53+0.04 
0.60+0.05 
0.28+0.03 
0.80+0.07 
0.44+0.04 
0.61+0.06 
0.67+0.06 
1.03+0.09 
1.36+0.12 



2.99+0.22 
13.14+0.95 
5.43+0.39 
27.67+1.96 
5.30+0.37 
8.77+0.63 
23.51+1.66 
29.68+2.10 
3.67+0.26 
6.27+0.44 
3.81+0.29 
25.19+1.78 
1.57+0.12 
16.27+1.17 
4.70+0.90 
12.00+2.40 
9.18+0.65 
9.40+0.69 
2.18+0.16 
2.25+0.17 
10.41+0.74 
3.21+0.23 
4.24+0.30 
2.12+0.16 
1.21+0.10 
1.90+0.14 
1.65+0.12 
1.72+0.13 
0.77+0.07 
2.21+0.16 
1.15+0.09 
1.25+0.12 
1.83+0.14 
2.54+0.19 
4.00+0.31 



7.05+0.50 
28.45+2.02 
13.38+0.94 
68.36+4.80 
12.56+0.88 
19.84+1.41 
56.24+3.95 
70.70+4.96 
9.31+0.66 
14.97+1.05 
8.71+0.63 
61.21+4.30 
3.40+0.25 
35.37+2.49 
10.80+2.20 
30.80+6.20 
22.97+1.62 
21.79+1.57 
5.31+0.38 
3.89+0.28 
21.92+1.54 
7.88+0.56 
10.54+0.74 
4.23+0.30 
2.58+0.19 
4.41+0.31 
4.20+0.30 
4.15+0.30 
1.64+0.12 
4.89+0.35 
2.55+0.19 
2.23+0.18 
4.06+0.29 
5.24+0.37 
8.92+0.68 



11.43+2.29 
34.80+6.96 
22.79+4.56 
117.3+23.5 
19.56+3.91 
25.81+5.16 
95.96+19.2 
94.58+18.9 
17.37+3.47 
24.10+4.82 
12.33+2.47 
94.14+18.8 
5.12+1.03 
41.18+8.24 
20.50+4.10 
51.10+10.2 
35.95+7.19 
30.80+6.16 
9.91 + 1.98 
4.08+0.82 
27.06+5.41 
12.78+2.56 
21.01+4.20 
6.87+1.38 
4.45+0.89 
6.34+1.27 
7.64+1.53 
6.13+1.23 
2.48+0.50 
6.12+1.22 
3.98+0.80 
2.03+0.41 
5.97+1.20 
8.12+1.63 
(...) 



10.10+2.02 
22.99+4.60 
19.85+3.97 
100.3+20.0 
13.98+2.80 
16.63+3.33 
90.52+18.1 
70.32+14.1 
18.46+3.69 
18.45+3.69 
10.65+2.13 
72.34+14.5 
4.13+0.83 
23.76+4.76 
16.10+3.20 
46.10+9.20 
28.35+5.67 
20.38+4.08 
9.72+1.95 
2.22+0.45 
15.16+3.04 
9.84+1.97 
21.76+4.35 
5.99+1.20 
4.95+0.99 
4.42+0.89 
6.36+1.27 
4.37+0.88 
1.74+0.35 
2.96+0.60 
2.81+0.56 
0.63+0.14 
3.90+0.78 
6.30+1.26 
(...) 



7.6+0.1 
7.9+0.1 
7.3+0.1 
8.0+0.1 
7.4+0.1 
7.7+0.1 
7.9+0.1 
8.1+0.1 
7.1+0.1 
7.4+0.1 
7.3+0.1 
8.0+0.1 
6.9+0.1 
8.0+0.1 
7.3+0.1 
7.6+0.1 
7.6+0.1 
7.6+0.1 
7.5+0.1 
7.3+0.1 
7.8+0.1 
7.1+0.1 
7.7+0.1 
7.1+0.1 
6.7+0.1 
7.2+0.1 
6.8+0.1 
6.9+0.1 
6.6+0.1 
7.4+0.1 
6.8+0.1 
7.4+0.1 
7.2+0.1 
7.1+0.1 
7.1+0.1 



22.0+0.9 
18.7+0.7 
22.4+0.9 
22.8+0.9 
21.0+0.8 
19.5+0.7 
22.4+1.0 
20.6+0.7 
23.4+1.2 
21.3+0.8 
20.5+0.9 
21.2+0.8 
20.2+0.9 
18.4+0.6 
21.8+1.0 
22.9+0.9 
21.6+0.7 
20.2+0.7 
22.8+1.2 
15.6+0.8 
18.2+0.6 
21.7+0.9 
23.7+1.3 
19.0+1.1 
20.9+1.5 
20.3+0.8 
23.2+1.0 
20.9+0.8 
19.9+0.9 
18.5+0.5 
20.0+0.9 
14.9+0.5 
19.8+0.8 
19.4+0.9 
22.7+3.6 



0.59 
0.68 
0.88 
1.36 
1.20 
1.28 
0.29 
2.28 
0.84 
0.86 
0.51 
1.09 
0.17 
1.13 
0.22 
0.12 
1.64 
1.42 
0.36 
2.73 
0.56 
1.34 
0.36 
2.79 
2.15 
1.14 
2.24 
1.69 
0.07 
2.57 
0.11 
0.52 
0.39 
0.87 
0.01 



The slope is closer to unity if we consider only galaxies of the 
M-sample, likely because this sample is more uniform in galaxy 
morphology (mostly Sc). The best correlation with the dust mass 
is found for the total gas mass, with a Pearson linear correlation 
coefficient of 0.91 for the M-sample and of 0.84 for the whole 
sample. In the last case the best-fit linear relation between M^,,,, 
and Mgas reads: 

logM,,, = 0.75(+0.09) X logMrf,,,, + 3.9(+0.6) (5) 

which implies a lower dust-to-gas ratio for galaxies with a lower 
gas content. If replaces in the computation of Mg„s the 
slope is 0.85+0.08 and r = 0.88. The decrease of the dust-to- 
gas ratio, as the gas mass decreases, is not due to the environ- 
ment: large galaxies, with a large gas mass, are not more HI de- 
ficient than small galaxies (with a small gas mass). In Figure |4] 
we show the relation between dust and gas mass for HI defi- 
cient galaxies with def(HI)> 0.5. The slope of the dust-to-total 
gas mass relation is 0.76+0.08, compatible with that found for 
the whole sample or for non-HI deficient galaxies (0.62+0.11). 
As expected, the dispersion in the M,/„5,-Mh/ relation decreases 
for HI deficient galaxies, which have dispersed the outer gas in 
the intercluster medium (bottom panel of Figure |4}. In the mass 
range covered by our sample, environmental effects on the dust- 
to-gas ratio (see next Section) do not depend on galaxy mass. 



Variations in the dust-to-gas ratio with the global gas content 
are instead due to internal evolution. Smaller galaxies in the S- 
sample have on average a higher HI content: they have been un- 
able to increase the dust abundance and to convert much of their 
atomic gas into molecular gas due also to their extended outer 
HI disks. As shown by Pappalardo et al. (2012) in a spatially re- 
solved study, the relation between dust emission and gas surface 
density in the inner regions is very tight, with a slope closer to 
one for metal rich galaxies. 

In the following we convert the flux density in an Herschel 
band, F^, into luminosity L,, = 47rD^vFv and express Ly in solar 
units. Hence we convert the CO J=l-0 total flux in Jy km s"' 
into Leo in solar units by multiplying Fco by 0.12 f^^^. The 
250 /im luminosity, similar to M,/,,,,, correlates better with the 
total gas mass rather than with the HI or H2 mass alone. 

The scatter in the 'Lco-^B,H,K,Ha relation, the B,H,K-band 
and Ho- luminosities, which we corrected for internal extinction, 
is larger than observed between Leo and the luminosities in the 
various Herschel bands. We display Leo versus L/za and M» in 
Figure |5] and below we give the relations between Leo and L^, 
Lh», M*, L250 for the whole sample (in solar units): 

log Leo = 1-10 X log L^ -6.7 r = 0.78 (6) 

logLeo = 0.87 X log L„„- 5.4 r = 0.83 (7) 
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Fig. 1. Fluxes in the 5 Herschel bands (in Jy units) versus the total CO J= 
M-sample, open circles galaxies in the S-sample. 



1-0 flux (in Jy km s '). Filled circles are galaxies in the 



log Leo = 0-78 X log M. - 4.3 r = 0.85 
log Leo = 1 -07 X log L250 - 5.6 r = 0.92 



(8) 
(9) 



The correlations of Mgas with Lho or M, are also less tight 
(see Figure |5]l than with L250 and much weaker correlations are 
found for Mhi- This confirms that the main correlation with the 
CO J=l-0 line intensity is indeed established by the cold dust 
emission. In our sample, which spans a limited metallicity range, 
the CO J=l-0 line emission traces molecular hydrogen which 
forms on dust, and this explains the correlation. Star formation 
is not ubiquitous in molecular clouds and, as a result, a less tight 
correlation of the CO J=l-0 line emission with Ha luminosity is 
observed. 



4. Environmental effects on the dust and gas 
content 

One aspect related to the environment which we are going to 
address is whether the global dust-to-gas ratio can be used 
as a metallicity indicator for cluster galaxies. We devote the 
rest of this section to investigate if the environment can en- 
hance or inhibit the persistence of dust and molecules in galax- 



ies (Star k et al.iri986HFumapani & Gavazzill2008h . The environ- 
mental efi^ects are examined primarily in terms of HI deficiency 
which we have comput ed for all galaxies in our sa mple follow- 
ing the prescription of iHavnes & Giovanellil(ll984l) . 

In Figure jS^a) we show the consistency of the metallicities 
of the central regions of galaxies in our sample with the mass- 
metallicity relation. In Figure^/?) the dust-to-gas ratio is plotted 
as a function of stellar mass. The dashed line shows the mass- 
metallicity relation scaled so that its peak reproduces the aver- 
age dust-to-gas ratio found for the M-sample, i.e. -1.9+0.15 in 
log. This implies a dust-to-gas ratio of 0.013, a facto r 1.8 higher 
than the value of 0.007 found for the Milky Way ( Draine et al.l 
2007), whose central regions have metal abundances comparable 
to those of the M-samp le, but only sligh tly higher that the theo- 
retical value of 0.01 (Draine et al.ll2007l) . In panel (c) we show 
the better agreement with the mass-metallicity relation if MJ^ is 
used instead (i.e. a constant CO-to-Ha conversion factor). In this 
case the discrepancy between the average dust-to-gas ratio in the 
Milky Way and that of the M-sample is reduced. 

Galaxies in the S-sample lie mostly below the mass- 
metallicity relation both in panel (b) and (c). For galaxies in 
the S-sample the dust-to-gas ratio is lower than what a lin- 
ear scaling with the metallicity would predict, similar to what 
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Fig. 2. Lower panel: The CO J=l-0 flux, in Jy km s measured 
in one beam at the galaxy center, versus the 250 /^m flux in Jy 
in the same area for galaxies in the M- and S-sample. The filled 
squares (red in the on-line version) are for beams of 100 arcsec 
FWHM, while the filled triangles (blue) are for beams of 45- 
55 arcsec FWHM. Upper panel: same quantities as in the lower 
panel but integrated over the whole galaxy. Filled circles are for 
galaxies mapped in CO while open circles are for galaxies whose 
total integrated CO flux has been inferred using pointed obser- 
vations (see text for details). 



Im unoz-Mateos et a find. These results holds even if we 

consider HI deficient and non-HI deficient galaxies separately. 
Likely, this is a consequence that metal abundances and dust- 
to-gas ratios sample different galaxy regions. The SDSS survey 
measures metal abundances at the galaxy center while the dust- 
to-gas ratio in this paper is a quantity integrated over the whole 
disk. Unenriched extended gaseous disks are more pronounced 
in low mass blue galaxies. The average dust-to-gas ratio for the 
M- and S-sample together is -2.1+0.3 in log. The best linear fit 
to the M- and S-sample data for Mj/Mga^ as a function of M, 
reads: 



logMrf„,,/M,,„, = 0.26 X logM. - 4.9 



(10) 

0.65) is shown by 



This relation (with correlation coefllcient r 
the dot-dashed line in Figure [SJZ?). 

To check whether the HI deficiency affects the dust-to-gas 
ratio we show this last variable as a function of the HI defi- 
ciency, def(HI), in the bottom panel of FigureQ We can see that 
as the HI deficiency increases, the dust-to-gas ratio increases but 
the correlation is weak (best fit slope is 0.26+0.06, r - 0.58). 
This trend can be due to two possible effects: the first one is 
gas removal which affects the dust content less because dust 
lies more deeply in the potential well of the galaxy. The sec- 
ond one is an evolutionary effect: if the galaxy entering the clus- 
ter turns off gas infall, it uses its gas reservoir to make more 
metals and stars. If this were the dominant effect, galaxies with 
higher metallicity (with respect to the metallicities predicted by 
the mass-metallicity relation) should be more HI deficient. We 



Fig. 3. The relation between dust and gaseous masses (HI,H2,and 
total gas mass in solar mass units) for the whole sample. Filled 
and open circles are for galaxies in the M-sample, and S-sample, 
respectively. The dashed line in the upper panel is the best-fit 
linear relation for the data shown. 
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Fig. 4. The relation between dust and gaseous masses (in so- 
lar mass units) for galaxies with def(HI)> 0.5. Filled and open 
circles are for galaxies in the M-sample, and S-sample, respec- 
tively. The dashed line in the upper panel is the best-fit linear 
relation for the data shown. 



do not find any correlation between the HI deficiency and the 
scatter in metallicity with respect to the mass-metallicity rela- 
tion. Hence, gas removal is responsible for the observed trend. 
The small panels on the right side of Figure |7] show the average 
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Fig. 5. The relations between the extinction corrected Ha lumi- 
nosity and stellar mass and the CO luminosity or total gas mass. 
Luminosities and masses are in solar luminosity and mass units, 
respectively. Symbols for galaxies in different samples are coded 
as in Figure [1] 




M. 



Fig. 6. In (a ) the mass-metallicity relation derived by 
iTremonti et al.l (l2004l) (dashed line) is plotted together with data 
for galaxies in our sample with a known SDSS metallicity, given 
in Table 1 . In panels (b) and (c) the mass-metallicity relation has 
been scaled so that its peak reproduces the average dust-to-gas 
ratio for the M-sample. In panel {b) is used to compute M^o, 
and the dot-dashed line is the best fit linear relation to galaxies 
in the M- and S-sample. In panel (c) M^^ is used in computing 
Mto,. 



values in 3 def(HI) bins. It is important to note that the dust-to- 
gas mass ratio does not increase as the galaxy becomes highly HI 
deficient, while between non-HI deficient galaxies and galaxies 
which have a mild deficiency this ratio increases by a factor two. 
This finding suggests that when the galaxy is highly disturbed 
the galaxy also becomes dust deficient. In terms of total mass, a 
small disturbance mostly affects the outer HI envelope and does 
not much affect the total dust mass since dust is more segregated 
in the inner regions of the disk (both in height and radially) . A 
similar conclusion has been reached bv lCortese et al.l (1201 2h by 
comparing a sample of field and cluster spirals. They find that 
the effect of dust removal from the star-forming disk of cluster 
galaxies is less strong than what is observed for the atomic gas. 

To investigate further whether some fraction of dust mass 
is lost as the galaxy moves through the cluster, we normalize 
the dust mass to the stellar mass. Dust removal is responsible 
for the correlation between the dust-to-stellar mass ratio and 
the HI deficiency shown in Figure [Tj/?). Galaxies with a higher 
HI deficiency have a lower dust-to-stellar mass ratio, only 25% 
of that observed in non Hl-deficient galaxies. The slope of the 
correlation is -0.48+0.08 and the Pearson linear correlation co- 
efficient is -0.72. Dust deficienc y in HI defici e nt V irgo spirals 
has been pointed out already bv ICortese et al.l (12010 ) in a spa- 
tially resolved study. The authors show that the extent of the 
dust disk is significantly reduced in HI truncated disk, due to 
stripping by the cluster environment There is the possibility that 
dust removal could be more effective in lo w-metalhcity, low- 
mass gal axies with shall ow potential wells (iBoselli & GavazzH 
[2006; B endo et alj|2007l) but the HI deficiency does not show 
any dependence on the stellar mass. 

To analyze possible environmental effects on the molecular 
component we use the Mn^-to-M, and M//,-to-Mgas ratios. We 
show these as a function of def(HI) in the top panels of Figure |7] 
The best fitting linear relation between log M//,/M, and def(HI) 
has a slope of -0.35+0. 11. The dependence on the HI deficiency 
of ratios involving the molecular mass, or the CO luminosity, 
is weaker (r~ 0.5) than for the dust component. If confirmed 
by future data this correlation implies that HI deficient galax- 
ies are also H2 deficient due to molecular gas removal or to 
quenching of molecule formation in perturbed gas. For our sam- 
ple the MH,/M,-def(HI) relation does not dependent on Xco- 
the slope of the correlation between the Lco-to-M» ratio and the 
HI deficiency is also shallow (-0.28+0.09) and consistent with 
that found for the M/f,-to-M» ratio. The gas molecular fraction, 
MnJMgas, increases as the HI deficiency increases (FigureQii)) 
because the atomic gas mass is more heavily removed in the 
cluster environment. 

As shown in Figure |7jc), highly HI deficient galaxies can 
quench or remove as much as 50-60% of their original molecu- 
lar content. If we consider only the mapped galaxies in the M- 
sample the decrease of H2 content is confirmed and becomes 
more extreme since the CO luminosity, as well as the H2 mass, 
shows a steeper decrease with the HI deficiency parameter A fu- 
ture comparison with a field galaxy sample is desirable to con- 
firm this finding. The cluster environment strips out part of the 
tenuous atomic ISM (e.g. outer disks and high latitude gas) more 
easily than dust and H2 gas, which is more confined to the mid- 
plane and to the inner star forming disk. In terms of total mass, 
dust and molecules are less affected by the environment than the 
HI gas. Highly disturbed galaxies have considerably lower dust 
content than their quiescent counterparts and higher molecular 
fractions. 

Figure [8] shows possible dependencies of galaxy properties 
on the distance from the cluster center Non-HI deficient galax- 
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Fig. 7. The dust-to-gas ratio, dust-to-stellar mass ratio, molecular-to-stellar mass ratio and the molecular fraction as a function of the 
HI deficiency. Panels to the left show the data with symbols as in previous figures: filled circles are galaxies in the M-sample, open 
circles are galaxies in the S-sample. Continuous lines in each panel are best fits to linear relation between def(HI) and the variable 
shown on the y-axis. Panels to the right of the figure show mean values in three HI deficiency bins. 



ies avoid the cluster center, where a higher fraction of Sab and 
Sb (open triangles) can be found. The dispersion across the 
whole cluster is large because we are using projected distances. 
The evidence that dust deficient galaxies and galaxies with a 
lower Lco/M, are located in projection near the cluster center 
is marginal (see bottom and top panel of Figure|8). 

5. Summary 

The relations between dust and gaseous masses have been ex- 
amined in this paper for a magnitude limited sample of late type 
Virgo Cluster galaxies in the HeViCS fields. The Herschel Space 
Observatory has detected cold dust emission for all galaxies in 
our sample with the PACS and SPIRE detector at 100, 160, 250, 
350 and 500 jum. We summarize below the main results: 

- There is a tight correlation between the total CO J=l-0 flux 
and the cold dust emission. Molecules are more closely re- 
lated to cold dust rather than to dust heated by star formation 
or to other galaxy integrated quantities such as blue luminos- 
ity, star formation or stellar mass. 



- Dust mass correlates better with the total 
(atomic+molecular) gas mass rather than with the atomic 
or molecular gas alone. Infrared-dim galaxies are unable 
to convert much of their atomic gas into molecular gas. 
Galaxies with a low molecular-to-total gas ratios have lower 
dust temperatures, are intrinsically faint and blue and a have 
a high gas-to-stellar mass ratios. 

- The dust-to-gas mass ratio typically decreases as the galaxy 
stellar mass decreases, faster than predicted by the mass- 
metallicity relation. We find a population of galaxies with a 
lower dust-to-gas ratio than what a linear scaling with metal- 
licity would imply. Likely this is because the total gas mass 
of a galaxy includes the unenriched extended disk. In fact 
these galaxies are HI rich with a less prominent stellar mass. 

- Galaxies with a mild Hl-deficiency have a higher dust-to-gas 
ratios than non Hl-deficient galaxies. This implies a higher 
efficiency of the cluster stripping process on the HI gas rather 
than on the dust, which is more confined inside the disk. 
However, the cluster environment is able to remove some 
dust from the galaxy disk as well. The dust mass per unit stel- 
lar mass decreases as the HI deficiency increases, and highly 
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Fig. 8. The dust-stellar mass ratio, the HI deficiency and the ra- 
tio between the CO luminosity and the stellar mass are shown as 
a function of the projected distance from the Virgo cluster cen- 
ter Galaxies are coded as in Figure [1] The open (red) triangles 
indicate galaxies of morphological type Sab or Sb. 

HI deficient galaxies can have up to 75% less dust than non 
HI deficient galaxies. 

- There is evidence of molecular hydrogen quenching or strip- 
ping in HI deficient galaxies. This evidence is found by con- 
sidering the molecular gas per unit stellar mass. If its am- 
plitude is confirmed by future data implies that molecules 
are missing in Virgo HI deficient spirals, but to a somewhat 
lesser extent than dust. The molecular fraction increases as 
the HI deficiency increases because the cluster environment 
affects molecules to a lesser extent than atomic gas. 

- Our results have been confirmed using various assumptions 
on the CO-to-H2 conversion factor We cannot constrain pos- 
sible dependencies of Xco from galaxy properties, such as 
the disk morphology or metallicity, due to the limited dy- 
namical range of these parameters in the sample examined. 
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